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ABSTRACT 

An attempt has been made to study an oscillatory MHD free and forced convective flow past a vertical porous 
plate in slip-flow regime with variable suction and constant heat flux. The temperature of the plate oscillates in time about 
a constant mean and a uniform magnetic field is assumed to be applied transversely to the direction of the flow. The 
magnetic Reynolds number is assumed to be small so that the induced magnetic field can be neglected. The solutions to the 
governing equations are derived by regular perturbation technique with Eckert number (E) as perturbation parameter. 

The expressions for the velocity field, temperature field, skin friction at the plate in the direction of the flow and 
the plate temperature are obtained in non-dimensional form. The amplitudes and phases of the fluctuating parts of the skin 
friction and rate of heat transfer (Nusselt Number) are obtained in non-dimensional form. The effects of the Hartmann 
number (M ), the frequency of oscillation! (0 ), suction parameter) A ) and the rarefaction parameter! h) on these fields are 
discussed and demonstrated with the help of graphs. 

KEYWORDS: MHD, Viscous, Oscillatory, Electrically Conducting, Incompressible, Heat Transfer, Slip-Flow, Eckert 
Number 

1. INTRODUCTION 

The problems of convective flows arising in fluids as a result of interaction of the force of gravity and density 
difference caused by simultaneous diffusion of thermal energy and chemical species have been studied by many authors 
due to applications of such problems in Geophysics and Engineering. In this regard, we may cite the works done by Bejan 
and Khair [1], Trivision and Bejan [2], Acharya etc al [3], Raptis and Kafousias [4], Ahmed etc al ( [5]) . Unsteady 
oscillatory free and forced convective flows plays a major rule in chemical engineering, turbo machinery and aerospace 
technology. The temperature change causes density variation which leads to the free convection in the fluid. The free 
convection flow is enhanced by superimposing oscillating temperature on the mean plate temperature. 

Soundalgekar[6,7,8] studied the free convection effects on mean velocity and temperature field of oscillatory 

flow past an infinite vertical porous plate with constant suction with or without transverse magnetic field. The analysis of 
viscous dissipation on the transient free convection flow past a semi-infinite vertical plate is studied by Soundalgekar 

et. al. [9] . The results of the combined free and forced convection flow of water at 4° C from a vertical plate with variable 

temperature was investigated analytically by Vighnesam and Soundalgekar [10] . Sahoo et. al [11] have investigated an 
unsteady MHD free convection flow of a viscous incompressible electrically conducting fluid past an infinite vertical 



www.iaset.us 



editor @iaset. us 



42 



J.K. Goswami 



porous plate subjected to constant suction and heat sink. The analysis of the MHD unsteady free convection flow past a 
vertical porous plate was made by Anwar [12] . Hussain et. al. [13] studied the effects of a fluctuating surface 
temperature and concentration on a natural convection flow from a vertical plate. 

Analytical solutions to the problems of the transient free convective viscous incompressible flow past a vertical 
plate with periodic temperature and variable suction in a slip-flow regime are presented by many authors. Some of them 

are Sharma and Choudhary [14] , Sharma[l5] , Jain and Sharma[l6] and Ahmed and Goswami[17]. This paper presents 

the extension work of Ahmed and Goswami [17] to investigate the effects of various parameters namely Hartmann 

number ( M ), the frequency of oscillation! CO ), suction parameter! A ) and the rarefaction parameter! h) on a flow with 
mass transfer in a slip-flow regime. 

2. MATHEMATICAL FORMULATION 

We now consider a two dimensional boundary layer flow of an oscillatory MHD free and forced convective flow 
past a vertical porous plate in slip-flow regime with variable suction in presence of a transverse magnetic field by making 
the assumptions. 

I. All the properties except the density in the buoyancy force term are constant. 

II. The Eckert number E is small. 

III. The Magnetic Reynolds number is so small that induced magnetic field can be neglected. 

IV. The plate is subjected to a normal suction velocity, which varies periodically with time about a constant mean 
v 0 (say). 

V. The Magnetic dissipation term in the energy equation is negligible. 

We introduce a coordinate system (x,y,z) with X-axis vertically upward along the plate, Y-axis perpendicular 
to it and directed in to the fluid region and Z-axis perpendicularly to XY plane. Let q = III + vj be the fluid velocity at 
the point ! ( A' , V , Z ) and B — B 0 j be the applied magnetic field with strength B lt , i and j being the unit vectors along 

X-axis and Y-axis respectively. Since the plate is of infinite length therefore all the physical quantities except the pressure 
p are independent of x. Under these assumptions the physical quantities are functions of y and t only. 



The equations governing the flow are 

Equation of Continuity 

I " 0 



=> v=-v 0 (l + £Ae imT ) 



( 2 . 1 ) 



Momentum equation: 
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dT 0 

Energy Equation 



a “-v 0 (l + eAe mT )^ = g/3(T-TJ + gj3(C-CJ+v^ + -B 0 2 (U-u) 
dy dy~ p 



( 2 . 2 ) 



c. 



dT 



sr.dT 



— ~v 0 (l + £Ae ) — 
dt dy 



k d 2 T 



2 - fdu^ 

+v 



pdy- 



(2.3) 



Species continuity equation 



dc n . imT dC _ d 2 C n d 2 T 
3F v ° il+eAe } 3y D «dy 2 T 3y> 



(2.4) 



The boundary conditions are 



u =h 



^ du ^ 



, F = T a +e(T-TJe u " , C = C+eiC-CJe' 6 ’' at y=0 



(2.5) 



u — > U , T — > F, , C — > at y — > ° o 

We introduce the following non-dimensional quantities 



( 2 . 6 ) 



yv 0 fv () 2 u _ T -T x 7/ C-C cja 

y = ^- 9 -, t = — — , co = — — , u = — , 0 = = — ==-, 7/ = — ,d = = — = 

^ ^ v 0 v 0 r^-71 v 0 C w -C. 



v 



D 



M 



VO 

, S 0 = — , v = — . 
0 A P 



vg/3(C w -CJ D _pC J , 
k 



G r = gfiv^^, G m = rs "^* -~\ P = !ZZJL, £ = _ 2° ^ = - 



C p <T m -TJ 



vji_ 

V 



PV o 



Where G r is the Grashof number for heat transfer, G m is the Grashof number for mass transfer, P is the Prandtl 
number, E is the Eckert number, h is the rarefaction parameter, k is the thermal conductivity, A is the suction 
parameter, V is the kinematics viscosity, M is the Hartmann number, C is the specific heat at constant temperature, f5 

is the coefficient of thermal expansion, S is the Schmidt number, CO is the frequency parameter and the other symbols 
have their usual meanings. 

The non-dimensional form of the equations (2.2), (2.3) and (2.4) are 



du 

~di 



-(1 + eAe‘“) 



du 

dy 



d 2 

= G r 6 + G m </> + — r + M(U —u) 
dy~ 



(2.7) 



30 

dt 



-(1 + eAe‘“) 



30 

dy 



1 d 2 0 ( du^i 

Pdy- ^3yJ 



2 . 8 ) 
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d 0 (l + eAe**)^- 1 d ^ + S dl ° 
S, " + £Ae ’dy S,dy ! "dy 1 


2.9) 


Subject to the boundary conditions 




u = h — , 0 = 1 + £e ica , (j) = \ + £e i0 * Xiy =0 
dy 


2.10) 


u = U , 0 = 0 , (/> = 0 at y — > 


2.11) 


3. METHOD OF SOLUTION 




Assuming the small amplitude of oscillation ( £ D 1 ), we represent the velocity U 
the plate as follows 


and the temperature 0 , near 


u(y,t)= u Q (y) + £e"°'u l (y) + 0(£ 2 ) 


(3.1) 


0(y,t) = 0 o (y) + £e ka 0 l (y) + o(£ 2 ) 


(3.2) 


</>(y, t ) = fo(y) + £e ,tt *fa(y) + o (£ 2 ) 


(3.3) 


2 

Substituting the above assumptions in (2.6) to (2.10) and equating the harmonic terms and neglecting £ , the 
equations for U Q , U l , 0 Q and 0 X with reduced boundary conditions are as follows. 


u 0 + u 0 -Mu 0 = —G r 0 o - GJ 0 - MU 


(3.4) 


«, + - (M + icd)u x = —G r 0 j - GJ X - Au 0 


(3.5) 


0 O + P0„ = -EPu () 2 


(3.6) 


0 x + P0 1 - iP(O0 x = -2PEu 0 «, -AP0 o 


(3.7) 


00 + S c0O = ~ S 0 S Ao 


(3.8) 


<Pl + S c ti - iPaxj) x = -AS c </> () - S o S c 0 x 


(3.9) 


Subject to the boundary conditions: 




Uq = h , M| = h , 0 Q = 1 , 0 X = 1 , </> 0 = 1 , (j) x = 1 at y = 0 

dy dy 


(3.10) 


u 0 = U , u x = 0 , 0 O = 0 , 0 X = 0 , 0 Q = 0 , </> x = 0 at >’ -A °o 


(3.11) 


To solve the equations we assume 




u 0 = a 01 + Eu 02 + 0(E 2 ) ; W, = u n + Eu l2 + 0(E 2 ) 


(3.12) 



Impact Factor (JCC): 2.0346 



NAAS Rating: 3.19 




Oscillatory MHD Free and Forced Convective Flow and Mass Transfr Through a 
Vertical Porous Plate in Slip-Flow Regime with Variable Suction and Constant Heat Flux 



45 



0 O — 0 ol + E0 Q1 + 0(E ) ; 0^ — 0 ^ j + E0^ + 0(E ) (3.13) 

0 o =0 oi +E 0 O2 + O(E 2 )-, $ l =$ ll +E$ l2 + 0(E~) (3.14) 



Substituting above in equations (3.4) to (3.11) and equating the terms independent of E and the coefficient of E 
in each equation and neglecting E~ , we obtained the following equations with respective boundary conditions are as 



follows. 

M oi + M 01 ~ Mu 0l = — 1 G r 0 0l — G m </> 0l — MU (3.15) 

m 02 m o2 — Mu 02 = ~G r 0 O2 - G m 0 O2 (3.16) 

u n + u n ~{M +iCD)u n = ~G r 0 n ~GJ) n -Au 0l (3.17) 

u l2 + u n — (M + id))u l2 = -G r 0 12 -G m <j) l2 -Au O2 (3.18) 

0 o ;+p0 o ;= o (3.i9) 

0 O2 +P0 O2 =~P(u ol ) 2 (3.20) 

0 n ' +P0 n -iO)P0 u = —AP0 OI (3.21) 

0 l2 + P0 l2 - iO)P0 u = -2 Pu m u n - AP0 O2 (3.22) 

K ^oi — — ^o^c^ot (3.23) 

^02 ^c^02 = —S()S C 0O2 (3.24) 

fy\\ +S Ai — iwSA i = ~AS c </)q X — S {) S c 0 u (3.25) 

*P\i *^c^i2 ~iO)S c (j) n = -AS c (j) Q2 — S o S c 0 l2 (3.26) 



Subject to the boundary conditions 



u 0l = h ^ M()l , u u =h ^-±L ; Um - fa ^ U ° 2 



dy 



dy 



dy 



l 12 



= h 



^^12 

dy 



^oi =1 ^02=0- =1.^12 =0. $>i=M)2=0, </> u =l, = 0 



at y = 0 



u Q1 U , u n 0, u Q2 0 ,u l2 O,0 m 0 ,0 Q2 d,0 n d,0 x2 0, </> w d,(f>Q 2 0, <j) n 0, 



(j) n = 0 at y — > oo 
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Solving the equations from (3.15) to (3.26) with subject to the boundary conditions we get 
0 M =e~ Py (3.27) 

0 U = B i e~ A "' + B 2 e~ Py (3.28) 

u M = \ x e~ Ky + \ 2 e~ Py + \ 3 e~ Scy + U (3.29) 

p -f p~ Py — f p~ 2 ^ y — f P ~ 2p y - f p~ 2Scy — f e~^ +P)y — F e~ {P+Sc)y — f p~^ +Sc)y 

w 02 1 01 e 1 02 c 1 03 e 1 04 e 1 05 c 1 06 c 1 07 c 

^=(1 + D l0 )e~ s ' y -D l0 e- Py (3.31) 

</> n = D n e~ Ly + D n e~ Scy + D l3 e~ Py + D u e~^' y (3.32) 



j Z> 02 = D is e~ s ‘ y + D l6 e~ Py + D l7 ec~ 2Ay + D n e~ 2Py + D 19 e~ 2S ‘ y + D 20 e~ (p+A)y + D 2l e~ (P+s ‘ )y + D 22 e~ (A+s ‘ )y 

(3.33) 



u — T + T p~ ,Jy +T p~ 1/ky +T p~ ziy +T p^ cy +T p 

“02 *00 c ~ J -01 C ^- t 02 e ^ ± 03 c ^-'04 c ^-'05 e 



o~ 2A i- v . 



- 2 S c y 



-(A +P)y 



+ T 06 e 



~(P+S c )y 



+ T 01 e 



-(\+S c )y . rp -S 

" r 1 08 e 

(3.34) 



«, , = V-*’ + T, + V"‘-’ + + V 



-Ay 



(3.35) 



0 n =H oo e~ A ' y +H n e HAi+h)y +H l2 e~ (Ai+A " )y +H l3 e~ (A ' +P)y +H u e~ (A ' +S)y + H 15 e~ iA+L)y 



+H,,e~ 2Aty + H.^e 



~(p+h)y 



+ H l& e 



~(P+\ i)y . 



+ H l9 e~ 2Py + H^e 



~(S e +P)y 



20 c 



+ H 2l e 



-(P+L)y 






~(s c +h)y 



+ H 23 e 



-(S c +Ai)y 



+ H 24 e- 2s ' y + H 25 e 



~(S C +L)y 



+ H 26 e 



-Py 



(3.36) 



^12 - A 3 e L ' + A 4 <? +D 25 e + D 26 <? 2/1 ' + £> 27 <? _P ' +Z) 28 e 









-2Ay 



,-2t>y . 



,-2S c y 



+ Z) 29 £ 



~(\+P)y 



+ D 30 e 



<S C +P)y 



+ D 2 ^€ 



-(A+A)y 



+ D-,^6 11 D^-,6 



-(\+k)y 



'32 



33 c 



+ D 34 e 



-(A+A i)y 



D 35 e~ (A+L)y + D 36 e~ (ll+ny + D 37 e~ (Ail+ny + D 3S e~ (P+L)y + D 39 e~^ +,1,y + D 40 e~^ +An,y + D 41 e 



,-(s c +h)y . 



-(S c +Ai)y 



~(S C +L)y 



(3.37) 



« 12 =V 00 e~' iy +V u e~™ +V„e~^ +,1,y + +V, A e~^ +ny +V 1K e 



,-Aiy 



-<A+/i)y I i/ -(A+Ai)y 

- 12<3 t v l 3 r 



(A+f)y , 1 / -(A+A)y 

14“ 



V l6 e~^ +L)y +V l7 e- 2Ay +V n e~ (p+h)y +V l9 e 



,-2Ay 



,-(i*+Ai)y 



+y 20 ^ 2/>) '+y 21 e- (P ^ )v +y 22 e- (P+i)y 



+y 23 e 



-(s c +4)y 



+ y 24 e 



-(s c +Ai)y 



+ V 25 e^ y +V 26 e-^ >y +V 27 e-" y +V 2S e^ y + V 29 e~*' y + V 30 e 



~(S C +L)y 



o~ p y 



-py 



•,-A.v 



-\y 



(3.38) 



Substituting equations (3.27) to (3.38) in (3.1) and (3.2) we get the expressions for velocity and temperature 
profiles as given below. 



u(y,t ) =u 0l (y) + Eu 02 (y) + e(u u + Eu l2 )(cos cot + i sin cot) 



(3.39) 



9{y, t) = 6 w (y) + E9 (n (y) + £{6 ^ , + Ed n )(c os cot + i sin cot) 



(3.40) 
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, 1 + Vl + 4M P + \Ip 2 +ia>P , 1 + Vl + 4M +iQ) 

Where /tj = , n x = , = 



R _i» R r - G r (l + Ph)-U(P 2 -P-M ) _ 

£>1 1 B 2 , P 2 l , B 3 /I . 7 i x/r.2 n ’ "4 n 2 



co ■ (1 + hA l )(P 2 -P-M) " P--P-M 

The other constant are obtained but not presented here for sake of brevity. 

4. COEFFICIENT OF SKIN-FRICTION 

The non-dimensional form of Skin-friction at the plate y = 0 is given by 



*o = 



du 

M Jy_ 

PJ 



^ du^ 



a >’A.=o 



/rv\ , _ i(j)t ' / r\\ _0 , _ ztMf _1 

= u 0 (0) + £e Mj (0) = r 0 + t 0 



y=o 



= A + £ I fi I cos(tw + a) 



, r {/, 

Where I B \- ^U R + U, , tan a = — 

Ur 

U R = Real parts of (t< n (0) + Eu n (0)^ , U 2 = Imaginary parts of (m u ( 0) + £w 12 (0)j 

The expressions for U R and U 1 are obtained, but not presented here for the sake of brevity. 

5. COEFFICIENT OF HEAT-TRANSFER 

The non-dimensional form of rate of heat transfer (Nusselt number N ) at the plate is given by 



k 


Ah 

ro 


1 


fde') 


i ' 


pv 0 C p {T m 


-TJ 








ii 

0 

1 

*0 

1 



N„ =■ 



= +£ I H I cos(cot + j3) 



^ i pc'™ 



dy 



y-0 



H, 



Where I H I = H R + H } and tan /3 = — — 

LL „ 



H R and H j are respective real and imaginary parts of 



d y J v-0 
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6. RESULTS AND DISCUSSIONS 

The graphs for amplitudelBI (versus co and h), the phase tana (versus co and h) of fluctuating part of the 
non-dimensional skin friction T and the amplitude IHI (versus co and h), the phase tanP( versus h) of the fluctuating part of 

the non-dimensional rate of heat transfer N have been presented. Throughout our investigation the Prandtl number P is 
taken to be equal to 0.71 which corresponds to the air, S=0.22 for /7 , , the free steam velocity U is equal to 1, the 
Grashof number for heat transfer G r = 5, the Grashof number for mass transfer G =2 and the Eckert number E is 
assumed to be 0.05 .The values of other parameters chosen arbitrary. 

Figures land 2 respectively demonstrate the behaviour of |fi against frequency (co) for different values of 

Hartmann number ( M ) and suction parameter ( A ). Both figures indicate that |/f sharply decreases as co increases, for 

small values of co and then it decreases steadily and slowly. Moreover an increase in the values of ( M ) causes I B I to 
decrease and the graphs are distinguished for small and moderate values of M, whereas opposite characteristics are noticed 
in the second figure. 

It is seen from figures 3 and 4 that the effects of Hartmann number ( M ), suction parameter (A) on 
I H I respectively are almost similar to the effects of these two parameter on I B I in figl and fig2. 

Figure 5 exhibits the variation of tana against h . It is observed that for CO < 5 the phase (tana) increases sharply 
and for large values of co it is insignificant It is also inferred from this figure that tana — > 0 as CO — > 00 irrespective of the 
values of M . 

Figures 6 and 7 are demonstrated the variation of tan a against h . It is observed from figure 6 that an increase in 
the values of h causes tan CX increases for small values of M and it remains constant for moderate values of M . It is 
also marked from figure 7 that tan CX decreases and the graphs are away from each other as h and M both increases. 
Further from figure 6 it is inferred that tan CX is distinguished initially for different values of M and as h — > 00 it is 
insignificant and these behaviour takes reverse trend in figure 7. 

Figure8 exhibits the variation of tanp under the influence of rarefaction parameter h and the Hartmann number 
M . It is clear from the figure that an increase in causes h tanP increases as M increases but for small values of M 
causes tanP remains constant. It is also noticed from the figure that there is no influence of rarefaction parameter h on tanP 
for small values o M . The same figure also shows that for a fixed value of h an increase in the values of M causes a 
growth in tanp. 
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Figure 1: Amplitude (IBI) of Skin Friction Versus Frequency (tl) When A =2 and h =0.4 




0.2 0.4 0.6 0.8 1 . 1.5 2 . 2.5 3 4 5 

frequency!®) 



Fig2: Amplitude of Skin-Friction IBI Against Frequency Parameter fi When M =1 




Figure 3: Amplitude of Heat Transfer (IHI) Against Frequency (fl) When A =2 and h =0.4 
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Figure 4: Fluctuation Part of Plate Temperature IHI Against Frequency Parameter f> When M =1 and h =0.4 
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Figure 7: Phase of Skin Friction (Tana) Against Rarefaction Parameter ( h ) When A =1 and C0= 3 




Figure 8: Heat Transfer Phase Tanp versus Rarefaction Parameter ( h ) When A =1 and <W=3 
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